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DYNAMICAL SYSTEMS ON LEIBNIZ 
ALGEBROIDS 



Gheorghe IVAN and Dumitru OPRI§ 



Abstract. In this paper we study the differential systems on Leibniz algebroids. 
We introduce a class of almost metriplectic manifolds as a special case of Leibniz 
manifolds. Also, the notion of almost metriplectic algebroid is introduced. 

These types of algebroids are used in the presentation of associated differential 
systems. We give some interesting examples of differential systems on algebroids 
and the orbits of the solutions of corresponding systems are established. 1 

1 Introduction 

Lie algebroids have been introduced repeatedly into differential geometry 
since the early 1950 's, and also into physics and algebra, under a wide variety 
of names. 

It is well known that there exists a one-to-one correspondence between 
Lie algebroid structures on a vector bundle 7r : E — > M and linear Poisson 
structures on the dual vector bundle ir* : E* — > M. This correspondence can 
be extended to much wider class of binary operations (brackets) on sections 
of 7r on one side, and linear 2 - contravariant tensor field on E* on the other 
side. It is not necessary for these operations to be antisymmetric or to satisfy 
the Jacoby identity. The vector bundle 7r togheter with a bracket operation, 
or the equivalent 2 - contravariant tensor field, will be called an algebroid. We 
mention the concept of Loday algebras, i.e. Leibniz algebras in the sense of 
Loday which are " non antisymmetric Lie algebras " . 

Weinstein 's paper on Lagrangian mechanics and groupoids (see [6]) roused 
new interest into the field of algebroids and groupoids. Weinstein introduces 
" Lagrangian systems " on a Lie algebroid by means of a Legendre-type map 
from E to E* associated to a given function L on E. The local coordinate 
expression of such equations reads 
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where the (x l ) are coordinates on M, (y a ) are fibre coordinates on E and the 
C^ b are structure functions coming from the Lie algebroid structure. Note that, 
more generally, the equations of the form 



were called " second - order equations on a Lie algebroid " by Weinstein. 

This paper shows that various relevant dynamical systems can be described 
as vector fields associated to smooth functions via a bracket that defines a 
Leibniz structure, defined of Leibniz algebroid. 

In this paper we go all the way in this direction and we work with a bracket, 
first introduced in (Grabowski and Urbanski, [2]), that is just required to be 
linear and a derivation on each its entres. The derivation property, also known 
as the Leibniz rule, justifies why we refer to this structure as Leibniz bracket. 

This construction should not be mistaken with the Leibniz structures (also 
called Loday algebras) introduced by Loday ([ 4 ]) in the algebraic context. 

The notion of Leibniz algebra is a certain noncommutative version of a Lie 
algebra. Leibniz algebroids are bundles with a local Leibniz algebra structure 
on its sections and a bundle map to the tangent bundle (called , anchor), which 
maps the bracket into the usual Lie bracket of vector fields. 

2 Leibniz systems. Almost Leibniz systems 

Let M be a smooth manifold of dimension n and let C°°(M) be the ring 
of smooth functions on it. A Leibniz bracket on M is a bilinear map [■,■] : 
C°°(M) x C°°(M) -> C°°(M) such that it is a derivation on each entry, that 
is, for all f,g,h e C°°(M) the following relations hold: 



x" = 



y a = f a (x,y) 



[fg, h] = [f, h]g + f[g, h] and [/, gh] = g[f, h] + [/, g]h. 



(1) 



We will say that the pair (M, [•,•]) is a Leibniz manifold. If the bracket [•, •] 
is antisymmetric, that is [f,g] = —[<?,/], (V)/, g G C°°(M), then we say that 
(M, [•,•]) is an almost Poisson manifold. 
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Let (M, [•,•]) be a Leibniz manifold and let h be a smooth function of M. 
From the fact that a Leibniz structure is a derivation there exists the vector 
field Xh on M uniquely characterized by the relation : 

X h (f) = [f,h] for any / G C°°(M). (2) 

We will call X h the Leibniz vector field associated to the Hamiltonian func- 
tion h G C°°(M). Note that since [•, ■] is a derivation on each their arguments, 
they only depend on the first derivatives of the functions and thus, we can de- 
fine the 2- contravariant tensor field B on M given by: 

B(df,dg) = [f,g], for all f,geC°°(M). (3) 

We call dynamical Leibniz system, the dynamical system associated to the 
vector field Xh- 

If (x l ),i = l,n is a system of local coordinates on M, in a local chart the 
dynamical Leibniz system is given by: 

• • dh 

x l = [x\ h], where [x\ h] = X l h = B{dx\ dh) = B 13 — . (4) 

Note that the symplectic and Poisson manifolds are particular cases of 
Leibniz manifolds. We now briefly present other non trivial examples. 

Example 1. {%) Pseudometric bracket and gradient dynamical sys- 
tems. Let g be a 2 - contravariant non degenerate symmetric tensor field on 
M. Let [•,•]: C°°(M) x C°°(M) -> C°°(M) be the Leibniz bracket defined by: 

[f,h]=s(df,dh), for all f,heC°°(M). (5) 

The bracket given by (5) is called the pseudometric bracket associated to g. 
This bracket is clearly symmetric and non degenerate and the Leibniz vector 
field Xh associated to any function h G C°°(M) is such that Xhf = g(df, dh). 
These brackets are also called Beltrami brackets. 

Let on R 3 the constant 2 - contravariant tensor field g = (g tj ) given by: 

/ S171 \ 
g = s 272 (6) 
V S 3 73 / 

where the parameters si, s 2 , s 3 G { — 1, 1} and 71,72, 73 are real numbers satis- 
fying the relation 71 + 72 + 73 = 0. This system happens to be a particular of 
pseudometric bracket and the differential system for h = x^x^x 2, is given by: 

i = S171X x , x — s 2 72^ x , x — S373X x . (7) 
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iii) The almost metriplectic systems . Let P be a 2 - contravariant 
antisymmetric tensor field on M and g be a 2 - contravariant non degenerate 
symmetric tensor field on M. We define the bracket [■, ■] : C°°(M) x C°°(M) — > 
C°°(M) by: 

[f,h] = P(df,dh)+g(df,dh), for all /,fter(M). (8) 

This bracket is clearly a Leibniz bracket. We will say that (M, P, g, [•, •]) 
is an almost metriplectic manifold. If the tensor field P is Poisson, then 
(M, P, g, [•, •]) is a metriplectic manifold, see [1]. 

The differential system associated to the bracket (8) is given by: 

dh dh 

x l = [x\ h], where [x\ h] = P{dx\ dh) + g(dx\ dh) = P^— + g lJ T—- (9) 

Let be the 2 - contravariant tensors fields P = (P* J ),g = (g^) on R 3 and 
the function h G C°°(R 3 ) given by: 

x 3 —x 2 
-x 3 x 1 \ (10) 
x 2 —x 1 



(— a\{x 2 ) 2 — a|(x 3 ) 2 aia 2 x 1 x 2 a^^x 3 \ 

a 1 a 2 x 1 x 2 — a 2 (x x ) 2 — a\(x 3 ) 2 a 2 a 3 x 2 x 3 (11) 

aia 3 x 1 x 3 a 2 a 3 x 2 x 3 — a 2 (x x ) 2 — a\{x 2 ) 2 J 

h = I( ai + ^(a; 1 ) 2 + 1(02 + l)(x 2 ) 2 + 1(03 + l)(x 3 ) 2 . (12) 

The differential system associated to the bracket (9) is given by: 

x 1 = (a 3 — a 2 )x 2 x 3 + a 2 (ai — a 2 )x l {x 2 ) 2 + a 3 (ai — a 3 )x 1 {x 3 ) 2 
x 2 = (ai — a 3 )x 1 x 3 + a 3 (a 2 — a 3 )a; 2 (x 3 ) 2 + a\{a 2 — ai)x 2 (a; 1 ) 2 (13) 
x 3 = (a 2 — ai)x 1 x 2 + ai(a 3 — ai)x 3 (x 1 ) 2 + a 2 (a 3 — a 2 )x 3 (x 2 ) 2 

The differential system (13) is the revised system of the rigid body , see [ 
1]. ' □ 

Let M be a smooth manifold and P and g be two 2 - contravariant tensors 
fields on M. We define the bracket [•,(•, •)] : C°°(M) x C°°{M) x C°°{M) — > 
C°°(M) by: 

[f,(h 1 ,h 2 )]= P(df,dh 1 )+g(df,dh 2 ), for all /,/ il ,/ i2 GC 00 (M). (14) 
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Proposition 1. The bracket (14) satisfy the following relations: 

[hf, (h, h 2 )\ = h[f\ {h u h 2 )] + /[/i, (h, h 2 )\ (15) 

[[/, Khi, h 2 )\ = l[f, (h u h 2 )] + h,P{df, oil) + h 2 g(df, dl 2 ) (16) 

forallf,f 1 ,h 1 ,h 2 ,leC°°(M). □ 
The bracket given by (14) is called the almost Leibniz bracket, and we say 

that (M, P, g, [•,(•, •)]) is an almost Leibniz manifold. 
By direct calculation we obtain the following result. 

Proposition 2. // the tensor field P is antisymmetric and the tensor 
field g is symmetric and there exist h±,h 2 G C°°(M) such that P(dh 2 ,df) = 
0, g(dh 1 ,df) = for all f e C°°(M), then the bracket (14) satisfies the 
relation: 

[f, (hi,h 2 )\ = If, (h, h)}, for h = h 1 + h 2 . (17) 

□ 

From the fact that the almost Leibniz structure (14) satisfies the relation 
(15) , there exists the vector field X hlfl2 on M uniquelly characterized by the 
relation : 

X hlh2 (f) = [f,(hi,h 2 )] for any feC°°(M). (18) 

We will call an almost Leibniz dynamical system, the dynamical system 
associated of the vector field Xh x h 2 given by (18). 

In system of local coordinates (x l ),i = l,n on M, the almost Leibniz 
dynamical system is given by: 

x l = [x\(hi,h 2 )} (19) 

where 

K(^M1=^ = ^|+^|. (20) 

Example 2. Let be the 2 - contravariant tensors fields P = (P lj ),g = (g lj ) 
on R 3 and the functions h x ,h 2 G C°°(R 3 ) given by: 

P= | -1 x 1 | , g= | -(x 3 ) 2 | (21) 
hi = \(x 2 ? + \(x*? , h 2 = \(x^ + x\ (22) 
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The almost Leibniz dynamical system associated to structure (R 3 , P, g, [•,(•, •)]) 
for the above tensors fields and functions is given by: 

•1 2-2 13-3 12/ 2\2 /oo\ 

x = x , x —xx, x = —x x — (x ) . (23) 

□ 

Example 3. If we take the tensors fields P = (P^) and g = (g lj ) on R 3 
and the functions h\,h 2 G C°°(R 3 ) given by: 

—x 3 x 2 

x 3 | , g = | j (24) 
-x 2 

h 1 = 1 -(xr + x\ h 2 = l -{x 2 ) 2 + \{x*) 2 (25) 

then the almost Leibniz dynamical system associated to (R 3 , P, g, [•, (•, •)]) for 
the above tensors fields and functions is given by: 

x = x , x = x x , x — —x x —xx. (26) 

□ 




3 Leibniz algebroids 

Let M be a smooth manifold of dimension n, let i\ : E — > M be a vector 
bundle and 7r* : E* — > M the dual vector bundle. By E = T(M, £") we denote 
the sections of 7r. 

A Leibniz algebroid structure (pseudo - Lie algebroid structure ) on a vector 
bundle it : E — > M is given by a bracket ( bilinear operation ) [•, •] on the 
space of sections £ and two vector bundle morphisms pi, p 2 : £" — > TM (called 
the Ze/it and the np/ii anchor, respectively) such that 

[fa 1 ,ga 2 \ = f pi{a 1 ){g)a 2 - gp 2 {(J 2 ){f)a 1 + fg[<Ti,a 2 ] (27) 

for all o-i,^ G £ and f,gE C°°(M). 

A vector bundle 7r : £" — > M endowed with a Leibniz algebroid struc- 
ture ([•,•], pi, P2) on E 1 , is called Leibniz algebroid over M and denoted by 

M, Pi, P2)- 

A Leibniz algebroid with an antisymmetric bracket [■, ■] ( in this case we 
have pi — p 2 ) is called pre - Lie algebroid. A Leibniz algebroid with a 
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symmetric bracket [•, •] ( in this case we have pi = —p\ ) is called symmetric 
algebroid. 

In the following, we establish a correspondence between Leibniz algebroid 
structures on the vector bundle n : E — > M and the 2- contravariant tensor 
fields on bundle manifold E* of the dual vector bundle n* : E* — ► M. 

For a given section a G E, we define a function i^.a on E* by the relation: 

i E *cr(a) =< cr(7r*(a)),a >, for a EE*, (28) 

where < •, • > is the canonical pairing between and 

Let A be a 2 - contravariant tensor field on E* and the bracket [•, -]a of 
functions given by the relation : 

[f,g] A = A(df,dg) (29) 

for all f,ge C°°(E*). 

For a given 2 - contravariant tensor field A on E*, we say that A is linear, if 
for each pair (fa, fa) of sections of n* , the function A(diE*fa, (He* fa) defined 
on E* is linear. 

Theorem 1. (Grabowski and Urbanski, 1997, [2]) For every Leibniz al- 
gebroid structure on n : E — > M with the bracket [•, •] and the anchors pi,p2 
, ^/iere ea;zsts an unique 2 - contravariant tensor field A on £"* swc/i f/iai the 
following relations hold: 

iE*[o-i,a 2 ] = [i E *o- u i E *CF2] a (30) 



= **/]a , 7T*(P2^)(/)) = [7r7,^^]A (31) 

/or a// (7, o-i, o- 2 G E and / G C°°(M). 

Conversely, every 2 - contravariant linear tensor field A on £"* defines a 
Leibniz algebroid on E by the relations (30) and (31). □ 

Let (x l ), i = l,n be a local coordinate system on M and let {e±, . . . , e m } 
be a basis of local sections of E ( dim M = n, dim E = m ) . We denote by 
{e 1 , . . . , e m } the dual basis of local sections of E* and (x l , y a ) ( resp., (x l ,£ a ) 
) the corresponding coordinate system on E ( resp., E* ). 

It is easy to see that every linear 2 - contravariant tensor field A on E* has 
the form: 

A = cU '-k 4 + pl i 18 is* - ^ k- (32) 
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where C% b , p\ a , p l 2a G C°°(M) are functions of (x l ). 

The correspondence between A and a Leibniz algebroid structure is given 
by the following relations: 

d d 
[e a ,e b ] = C d ah e d , pi(e a ) = p\ a — , p 2 (e a ) = p\ a -Q~- ( 33 ) 

We call a dynamical system on Leibniz algebroid ir : E — > M , the dynamical 
system associated to vector field X h with h G C°°(M) given by: 

= A(df, dh), for all / G C°°(M). (34) 

In a system of local coordinates, the dynamical system (34) is given by: 

L=[L,h] A , x l = [x\h} A (35) 

where 

f)h nh dh 

[Uh]K = C d M^ b + P \ a ^- , [x\h] A = -pl a — . (36) 

If the Leibniz algebroid is a pre - Lie algebroid, then the dynamical system 
(35) and (36) is given by: 

^ dh dh dh ^ ^ 



If the Leibniz algebroid is a symmetric algebroid, then the dynamical sys- 
tem (35) and (36) is given by: 

^ dh dh dh ^ ^ 

Example 4. Let the vector bundle 7r : E = R 3 x R 3 — > R 3 and tt* : 
= R 3 x (R 3 )* — > R 3 the dual vector bundle. We consider on £"* the 2 - 

contravariant linear tensor field A, the anchors p±,p2 and the function h given 

by: 

-£s£ 3 &x 2 
P = | 6x 3 -^x 1 | , pi = | x 3 | (39) 
-6^ 2 6a: 1 

0-10 

p 2 =[ 1 -x 1 I , and h(x,£) =x 2 & + x 3 &. (40) 
x 1 
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The dynamical system (35) and (36) associated to Leibniz algebroid (R 3 x 
R 3 , P, g, pi, p 2 ) for h = x 2 £ 2 + ic 3 6 is given by: 

(J: 

I 6 = 

The dynamical system (41) is the dynamical system associated to the Maxwel- 
Bloch system. The orbits of the system (41) represented in the coordinate 
systems Ox l x 2 x 3 and 0^ 2 ^ are given in the figures Fig.l and Fig. 2. 



x 6 (x 2 



-X 3 X 1 £ 1 

x 1 x 2 £ >1 



X = 

x 2 = 

X 3 = 



X 

x 1 x 3 
—x l x 2 



(41) 




Fig. 1 Fig. 2 



4 Almost metriplectic algebroids 

Let 7r : E — > M be a vector bundle , an antisymmetric bracket [•, -]i on E 
and [-,-]2 a symmetric bracket on E. A Leibniz algebroid endowed with the 
brackets [•, -]i and [•, is called almost metriplectic algebroid. 

Let Ai and A 2 the 2 - contravariant tensors fields on the bundle manifold 
E* associated with brackets [•, -]i and [•, -] 2 . 

We consider the bracket [•,(•,•)] : C°°(E*) x (C°°(E*) x C°°(E*)) -> 
C°°(£;*) given by: 

[f,{h 1 ,h 2 )]=A 1 (df,dh 1 ) + A 2 (df,dh 2 ), for all /, h±,h 2 G C°°(E*). (42) 

It is not hard to prove the following proposition. 
Proposition 3. The bracket given by (42) satisfy the relations: 



[fif, (hi, h 2 )\ = h[f, (huh*)] + /[/i, (fci, fca)] 



(43) 
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[/, l(h ± , h 2 )\ = l[f, (hi, h 2 )\ + h 1 A 1 (df, dl) + h 2 A 2 {df, dl) (44) 

for all f,fi,h,h 2 ,l e C°°(E*). □ 
From (43) it follows that X hlh2 : C°°(E*) x C°°{E*) -> C°°(£*) given by: 

= [/,(&!, for all /eC°°(r) (45) 

is a vector field on E 1 * associated to the functions hi, h 2 G C°°(E*). 

We will call an almost metriplectic dynamical system, the dynamical system 
associated to the vector field X hlh2 given by (45) with the bracket defined by 
(42). 

In a system of local coordinates on E* , the almost metriplectic dynamical 
system (45) is given by: 

L = X hlha (S a ) , x l = X hlh2 {x*). (46) 

where 

X hlh2 {£ a ) = Ai(d£ ,d/ii)+A 2 (df ,d/i2) , X hlh2 (x l ) = A 1 (dx\dh 1 )+A 2 (dx\dh 2 ) 

(47) 

Then the almost metriplectic dynamical system (46) and (47) is given by: 

t r<& C 9 hi , fid £ dh 2 _i_ n i dhi , A dh 2 

Sa — ^labSda^ ^2o6Sd ^ Pla ftr* " r > J 1a dx i ( ao\ 

■i _ n i dhi I n i dh 2 \^°) 

Proposition 4. Let A x and A 2 fre two 2 - contravariant tensors fields on 
E* for an almost metriplectic algebroid. 

(i) If there exist hi, h 2 G C°°(E*) such that A 1 (dh 2 , df) = and A 2 (dhi,df) = 
for all f G C ioo (E*), then for h = h + h 2 we have [f, (h, h)} = [f, (hi, h 2 )\. 

(ii) If the function hi is linear on the fibre of the vector bundle it* : E* — > 
M , then there exists A 2 with property that it is determined by hi. 

Proof, (i) From hypothesis A\(dh 2 ,df) = and A 2 (dhi,df) = for all 
/ G C°°(E*) and from (42) it follows that 

[/, (h, h)\ = Ai(df, dh) + A 2 (df, dh) = Ai(df, dhi) + M<V, dh 2 ) = [f, {hi, h 2 )\. 
(ii) From A 2 (df, dh) = for all / G C°°(E*) , we obtain: 
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If = iaK{x) from (49) follows 

dhj 



(50) 



Then, it is easy to verify that 

P % 2a = if a 7^ h P*2i = JJ ^2 k 



rid 
^2ab 



h i . 

. 1 i d]A 

h b l J 2a dx* 



dx k ' 



if a = i 
if a = 6 



□ 



satisfy the relations (50). 

Example 5. The dynamical system of the rigid body on the Leibniz alge 
broid tt : R 3 x R 3 — > R 3 . 

The dynamical system of the rigid body on R 3 is given by: 



x — p 



dh 



i,j = 1,2,3 



where h = |(ai(x 1 ) 2 + a 2 (x 2 ) 2 + a 3 (x 3 ) 2 ), ai > a 2 > a 3 > and 



P = (p ij ) 



—x 



x 




x 2 —x 1 



X 






Let 7r* : E* = R 3 x (R 3 )* -> R 3 the dual vector bundle of n : £ -> R 3 and 
(x l ,^j) be a coordinate system on £"*. Let the 2- contravariant tensor field Ai 
on £"* given by (Ci, — pi, pi) where pi — p and 

/ 

d = -e 3 x 3 o 

The dynamical system for Ai and hi(x,£) = aix 1 ^ + a 2 x 2 £ 2 + fl3^ 3 ^3 is 
given by: 



(51) 



The orbits of the dynamical system (51) represented in the coordinate 
systems Ox 1 x 2 x 3 and 0£i£ 2 £ 3 are given in the figures Fig. 3 and Fig. 4. 



f ?1 


= a 2 £ 3 x 2 £ 3 


- a^ 2 x 2 x 3 - a 2 £ 2 x 3 + a 3 £ 3 £ 2 




= a^ix 1 ;*; 3 


— ai^x 1 x 3 — a^x 1 + ai^ix 3 


< * 


= ai^^x 2 


- a 2 £ix 1 x 2 - ai^x 2 + a^x 1 






(a 3 - a 2 ):r 2 :r 3 


x 2 




(ai — as)^ 1 ^ 3 


( X 3 




(a 2 — ai)x x x 2 
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x 10~ 3 

>1(T 




Fig. 3 Fig. 4 



The tensor field Ai satisfies the relation Ai(o7i 2 , df) = for all / G C°°(E*) 
and h 2 (x, £) = x 1 ^ + x 2 £ 2 + x 3 &. 

Applying the Proposition 4 , it follows that there exists the 2 - contravariant 
tensor field A 2 with property A 2 (dh 1 , df) = 0, for all / G C°°(E*). 

The 2 - contravariant tensor field A 2 on E* is given by: (C 2 , p 2 , p 2 ), where 

/ - w l \ 

C 2 =\ %b-W 2 and 
V f £3 - V^s / 

(— Vi a\a 2 x x x 2 a\a^,x l x 3 \ 
aia 2 x 1 x 2 -V2 a 2 a 3 x 2 £ 3 
aia^x 3 a 2 a 3 x 2 x 3 — V 3 / 

with 

VI = (a 2 ) V) 2 + (a 3 ) 2 (x 3 ) 2 , ^ = (a 2 ) 2 x 2 £ 2 + (a 3 ) 2 x 3 £ 3 , 
V 2 = (aOV) 2 + (a 3 ) 2 (x 3 ) 2 , W 2 = K) + (a 3 ) 2 x 3 £ 3 , 
V 3 = (ax) V) 2 + (a 2 ) V) 2 , W 3 = (ai) + (a 2 ) 2 x 2 £ 2 - 
The dynamical system (46) and (47) is the following: 



f 6 


= (02(01 - 


- a 2 )x 1 x 2 — a^x 2 x 3 — a 2 x 3 )^ 2 + (a 3 (ai — a 3 )x x 


x 3 + a 2 x 2 x 3 + a 3 x 2 )£ 3 




= (oi(a 2 - 


- ai)x 1 x 2 + a 3 x 1 x 3 + ciix 3 )^ + (a 3 (a 2 — a 3 )x 2 


x 3 — a^x 3 — a^x 1 )^ 


< ^ 


= (02(03 


- a 2 )x 2 x 3 + aix 1 ^ 2 + a 2 x l )^ 2 + (ai(a 3 — a^x 1 
= (o 3 - a 2 )x 2 x 3 + a 2 (ai - a 2 )x 1 (x 2 ) 2 + a 3 (ai 


x 3 — a^x 2 — aix 2 )£i 




x 1 


- a 3 )x 1 (x 3 ) 2 




x 2 


= (ai - a 3 )x x x 3 + a 3 (a 2 - a 3 )x 2 (x 3 ) 2 + oi(a 2 
= (a 2 - ai)x x x 2 + ai(a 3 - a i )x 3 (x 1 ) 2 + a 2 (a 3 


-a 1 )x 2 (x 1 ) 2 
-a 2 )x 3 {x 2 ) 2 




X 3 



(52) 
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The orbits of the dynamical system (52) represented in the coordinate 
systems Ox l x 2 x 3 and 0^2^,3 for a± = 0.6, a<i = 0.4, a% = 0.2 are given in the 
figures Fig. 5 and Fig. 6. 




Fig. 5 Fig. 6 



Conclusion. In this paper are presented differential systems on Leibniz 
algebroids. The analysis of these systems realizes by studying the behavior 
of its solutions in the neighborhood of equilibria points. This analysis will be 
enterprised in the future papers. 
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